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ABSTRACT 

Models of the dynamical evolution of the early solar system following the dispersal of the gaseous 
protoplanetary disk have been widely successful in reconstructing the current orbital configuration 
of the giant planets. Statistically, some of the most successful dynamical evolution simulations have 
initially included a hypothetical fifth giant planet, of ice giant mass, which gets ejected by a gas giant 
during the early solar system’s proposed instability phase. We investigate the likelihood of an ice 
giant ejection event by either Jupiter or Saturn through constraints imposed by the current orbits 
of their wide-separation regular satellites Callisto and lapetus respectively. We show that planetary 
encounters that are sufficient to eject an ice giant, often provide excessive perturbations to the orbits 
of Callisto and lapetus making it difficult to reconcile a planet ejection event with the current orbit of 
either satellite. Quantitatively, we compute the likelihood of reconciling a regular Jovian satellite orbit 
with the current orbit of Callisto following an ice giant ejection by Jupiter of ~ 42% and conclude that 
such a large likelihood supports the hypothesis of a fifth giant planet’s existence. A similar calculation 
for lapetus reveals that it is much more difficult for Saturn to have ejected an ice giant and reconcile 
a Kronian satellite orbit with that of lapetus (likelihood ^ 1%), although uncertainties regarding the 
formation of lapetus, on its unusual orbit, complicates the interpretation of this result. 

Subject headings: methods: numerical — planets and satellites: dynamical evolution and stability 


1. INTRODUCTION 

Various solar system formation models argue that 
the giant planets underwent planetesimal driven mi¬ 
gration (e.g. Fernandez & Ip 1984; Malhotra 1995; 
Hahn & Malhotra 1999; Tsiganis et al. 2005) at early 
times (< 1 Gyr) following a dynamical instability. The 
Nice model, originally presented by Gomes et al. (2005), 
Morbidelli et al. (2005), and Tsiganis et al. (2005), with 
subsequent variants under the same name, has been 
the most successful in reproducing the settling of 
the four giant planets into their present orbital con¬ 
figuration (Tsiganis et al. 2005; Morbidelli et al. 2007; 
Levison et al. 2011), the Late Heavy Bombardment at 
^ 700 Myr (Gomes et al. 2005), the capture of Jupiter’s 
Trojan asteroids (Morbidelli et al. 2005), the capture of 
gas giant irregular satellites (Nesvorny et al. 2007), as 
well as the structure of the Kuiper belt (Levison et al. 
2008) and how its dynamical evolution led to the con¬ 
tamination of the outer asteroid belt by primitive trans- 
Neptunian objects (Levison et al. 2009). 

The precise nature of giant planet migration in the 
early solar system remains uncertain due to our lack 
of knowledge regarding each body’s initial conditions 
following their formation out of the solar nebula and 
the chaotic nature of the migration process. However, 
Morbidelli et al. (2009) argued that smooth divergent 
migration of the gas giants is unable to sufficiently excite 
their orbital eccentricities and inclinations to their ob¬ 
served values. Additionally, Brasser et al. (2009) showed 
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that such migration from an initial resonant conhgura- 
tion following the dispersal of the gaseous disk leads to 
excessive orbital eccentricities in the previously formed 
terrestrial bodies via sweeping secular resonances. A pro¬ 
posed solution, known as the jumping-Jupiter scenario 
(Brasser et al. 2009), invokes close encounters between 
the gas giants and an ice giant (IG) resulting in the step¬ 
wise migration of Jupiter and Saturn from their initial 
mean-motion resonance. This can sufficiently excite gi¬ 
ant planet eccentricities and inclinations whilst jumping 
over the problematic secular frequencies of the terres¬ 
trial planets. In addition, the jumping-Jupiter scenario 
does not disrupt the asteroid belt’s observed morphology 
(Morbidelli et al. 2010). 

A statistical study by Nesvorny (2011) of the dynami¬ 
cal evolution of the solar system during such a phase of 
frequent planetary encounters showed that the likelihood 
of reconstructing the current orbital conhguration of the 
four giant planets is increased when a fifth giant planet 
of approximately Uranian mass is included in the early 
solar system. The instability, which gives rise to mul¬ 
tiple planetary encounters, results in the ejection of the 
hypothetical fifth giant planet, reconstructing the outer 
solar system whilst preserving the orbits of the inner ter¬ 
restrial bodies over long timescales (Batygin et al. 2012; 
Nesvorny & Morbidelli 2012). Such planet scattering 
events (Rasio & Ford 1996; Weidenschilling & Marzari 
1996), applicable to any multi-body system, provide a 
potential explanation for the existence of the recently 
detected “free-floating” planets (e.g. Delorme et al. 2012; 
Liu et al. 2013; Luhman & Esplin 2014). 

In addition to reconstructing the current orbital con¬ 
figuration of both large and small bodies in the solar 
system, models attempting to achieve a full description 
of the solar system’s early dynamical evolution must re¬ 
quire the survival of the giant planets’ regular satellites 
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(e.g. Deienno et al. 2014). Regular satellites which are 
thought to form via accretion processes in circumplan- 
etary disks (Canup & Ward 2002; Mosqueira & Estrada 
2003), are expected to form on prograde, low-eccentricity 
orbits that are nearly coplanar with the host’s equato¬ 
rial plane and have relatively small semimajor axes. The 
current deviations of Callisto and lapetus’ orbits from 
circular, uninclined orbits therefore limit how close an 
IG could have come to the gas giants in the early solar 
system. 

In general, outer satellites, which are less tightly bound 
to the planet, will suffer larger perturbations during a 
close approach with an IG. In addition, the orbital eccen¬ 
tricities of these outer moons are only marginally damped 
through tides (which could otherwise mask the effects of 
early encounters). It is therefore the outermost regular 
satellites (Callisto around Jupiter and lapetus around 
Saturn) which provide the most stringent constraints. 

Specifically, Deienno et al. (2014) investigated whether 
the close encounters in the particular simulations of 
Nesvorny & Morbidelli (2012) (referred to as NM12), 
that best reproduced the giant planets orbital architec¬ 
ture, would excessively excite Callisto’s orbit. Given the 
interest in additional planets in the early solar system, 
we generalize this question to ask how likely it is to re¬ 
tain Callisto (lapetus) at its observed orbit following an 
ejection of an IG by Jupiter (Saturn)? Although in the 
jumping-Jupiter scenario only Jupiter may be responsi¬ 
ble for ejecting an IG, we include an analysis of close IG 
encounters with Saturn as we are more generally inter¬ 
ested in early solar system instability models in which 
either gas giant could undergo close encounters with the 
solar system’s IGs. It also permits the direct comparison 
of the likelihood of retaining a Gallisto-like satellite or¬ 
bit to the likelihood of retaining an lapetus-like satellite 
orbit following the ejection of an IG (see Sect. 6). 

In Sect. 2 we discuss the relevant properties of the 
satellites of interest. Sect. 3 summarizes our methods of 
investigation. Sect. 4 & 5 present our study’s results, and 
Sect. 6 presents our calculation of the likelihood of recon¬ 
ciling satellite orbits following planetary encounters with 
present-day orbits. We conclude with a detailed discus¬ 
sion in Sect. 7 and a summary in Sect. 8. 

2. SUMMARY OF SATELLITES: CALLISTO & lAPETUS 

Gallisto is the outermost Galilean satellite, moving on 
a nearly circular (e = 0.007) and uninclined {i ^ 0.28°) 
orbit. Callisto’s orbital period is ^ 16.7 days and is 
the only Galilean satellite not locked in a mean-motion 
resonance (Musotto et al. 2002). lapetus’ orbit is some¬ 
what more eccentric (e ~ 0.03) and circles Saturn every 
~ 79 days. Guriously, lapetus exhibits a significantly 
inclined orbit possibly due to inclined planetary encoun¬ 
ters between Saturn and an IG like those expected in the 
jumping-Jupiter scenario (Nesvorny et al. 2014). Gur- 
rent satellite orbital elements are summarized in Table 
1. These data were obtained from JPL HORIZONS.^ 

While Callisto does not participate in a MMR like 
the other Galilean satellites, its eccentricity evolution 
is nevertheless secularly coupled to that of the inner 
moons (Greenberg & Van Laerhoven 2011). The eccen¬ 
tricities of the inner satellites are more easily damped 

^ http://ssd.jpl.nasa.gov/?horizons 


TABLE 1 

Summary of Giant Planets and Satellites. 


Planet 


Jupiter 

Saturn 


Mp [Mq) 
Up (AU) 

Rp (AU) 

9.54 X 10-"^ 
5.20 

4.67 X 10-"‘ 

2.86 X lO-"' 
9.54 

3.89 X 10-4 

Satellite 


Callisto 

lapetus 


Ps (days) 

Q'S (-^p) 

es 

is n 

16.69 

26.93 

7.0 X 10-3 
0.28 

79.33 

61.15 

2.83 X 10-2 
7.49 


Note. — Data are from JPL HORIZONS: 

http: //ssd.jpl.nasa.gov/?horizons 

by tides than Callisto because of their smaller semima¬ 
jor axes and as a result of the secular coupling, ecaiUsto 
is damped on timescales faster than expected for an iso¬ 
lated planet-satellite system that is tidally locked. We 
take this into account as described in Sect. 4.2 follow¬ 
ing numerical simulations of this effect by Deienno et al. 
(2014). Conversely, the eccentricity damping of lapetus 
since the solar system’s instability phase, has been neg¬ 
ligible (Castillo-Rogez et al. 2007). 

3. METHODS 

We model the evolution of an IG heading for a close ap¬ 
proach with one of Jupiter or Saturn, and investigate the 
encounter’s effect on the gas giant’s satellites. We first 
determine viable encounter parameters that lead to the 
ejection of the IG. Within this set, we investigate if there 
are any satellites remaining around Jupiter (Saturn) with 
an orbit consistent with that of Gallisto (lapetus). 

3.1. Numerical Model: Simulating Planetary 
Encounters 

We consider a reduced-body subset of the solar system 
including the Sun, a satellite-hosting gas giant planet 
and an IG. The gas giant of mass Mp is initialized on a 
osculating, circular orbit with semimajor axis Op. By ap¬ 
proximating the gas giant’s orbit as circular, we remove 
any dependence on its orbital phase at the epoch of en¬ 
counter; tenc- Because the mass of the ejected IG is not 
well constrained, we select a fiducial value approximately 
equal to the mass of Uranus (M/ = 5 x as was 

used by NM12. A system of Keplarian satellites is placed 
in orbit around the gas giant (see Sect. 3.2 for a detailed 
description of satellite orbits). 

To limit the number of encounter parameters and com¬ 
putational cost of our survey, we assume a coplanar ge¬ 
ometry (we present a more detailed discussion of the ef¬ 
fect of inclined encounters in Sect. 7.1). Thus, at closest 
approach the IG’s velocity is perpendicular to its separa¬ 
tion vector from the gas giant and the encounter is fully 
specified at this time {tenc) by the impact parameter b, 
the relative planet velocity Vreh and the phase angle 9 
(see Fig. 1). 

We linearly step the impact parameter outward from 
bmin = 0.02 AU until encounters no longer lead to ejec¬ 
tions. Deienno et al. (2014) found that close encounters 
with b < bmin excessively excite the Galilean satellites; 
our estimates in Sect. 6 for the likelihood of retaining the 
observed satellite orbits following an ejection therefore 
represent a conservative overestimate due to our com¬ 
paratively gentle ejections. However, we note that values 







Could Jupiter or Saturn have ejected a fifth giant planet? 


3 



y 



Fig. 1.— Pictorial representations of a close planetary encounter between a gas giant (larger, red planet) and an unspecified ice giant 
(IG; smaller, blue planet) at the epoch of encounter tenc- The coplanarity of planetary encounters permits analysis solely in the xy-plane 
(z = 0). Left: the heliocentric reference frame with the Sun, at a distance of the gas giant’s semimajor axis, aligned with the x axis. The 
planet’s instantaneous velocity vectors are shown along with their trajectories (dashed lines) in the vicinity of tenc- Right: the reference 
frame is centered on the gas giant and is zoomed-in to depict the phase angle 9] the angle between the x axis and the IG’s position vector 
(vrel) tenc, measured in the counter-clockwise direction. Here, 6 ^ 310°. At tenc, the magnitude of the IG’s relative position vector is 
equal to the impact parameter b and is orthogonal to the relative velocity vector Vrel- In both diagrams, the dotted grey rings represent 
the ring of regular satellites in orbit around the gas giant. The scale used here is approximate as these schematics are not intended to be 
exact, but instead, are included to aid in the reader’s visualization of the experimental setup. 


smaller than bmin are less likely due to the reduced en¬ 
counter cross-section at small radii. We estimate the size 
of this effect by including a reduced simulation sample 
with b < bmin and find that our final results (Sect. 6) 
changed (fractionally) by < 12%. 

To determine which encounter parameters are capable 
of ejecting the IG, we uniformly sample 20,000 parameter 
combinations within b G [6mm,0.1 AU], Vrei G [l,5]uesc, 
and 9 G [0,27r), where Vgec = y/2GMp/b. We then 
remove any unphysical encounters in which the IG is 
unbound from the Sun prior to the encounter. In the 
cases of Jupiter and Saturn we find Ngim = 278 and 
Nsim = 274 valid parameter combinations respectively. 

To simulate an encounter with the aforementioned pa¬ 
rameters from Fig. 1, we integrate backwards in time 
until the absolute separation of the planets |r/ — Tpl > 2 
AU. At this separation a satellite with a Gallisto-like or 
lapetus like orbit feels a force from its host which is > 10^ 
times greater than that felt from the IG. The positions 
and velocities of all three massive bodies are then used 
as initial conditions for the forward simulations which 
include the satellites. These simulations are integrated 
forward in time towards the encounter at t = tenc and 
are halted at f = 2tenc- After 2tenc the influence of the 
IG on the satellites is again negligible and satellite orbits 
are no longer perturbed by the IG’s influence. In order 
for close encounters with 6 < O.I AU to be possible, nu¬ 
merous ‘soft’ encounters between the gas giants and IG, 
prior to tenc, were needed in order to build-up the IG’s 
eccentricity. These encounters will supply a small per¬ 
turbation to regular satellite orbits that, on average, in¬ 
crease satellite eccentricity over time (see Nesvorny et al. 
(2014) Fig. 4). However, the effect of encounters on 
satellite orbits is strongly dependent on 6, which is much 
larger than O.I AU for ‘soft’ encounters. We are there¬ 


fore only concerned with the strongest (final) encounter 
which leads to ejection of the IG as its effect dominates 
the final satellite orbits. 

Finally, we do not include the perturbations from the 
satellite host’s oblateness. This should be a good ap¬ 
proximation since the encounter timescale with the IG is 
much shorter than the precession timescales of the satel¬ 
lites due to the non-spherical shape of its host planet. 

3.2. Numerical Model: Initializing Satellites 

In each encounter simulation, we include a ring of Gal- 
listo or lapetus analog satellites around the gas giant 
planet. The satellite ring consists of TV = 100 non¬ 
interacting test particles with azimuthal positions ran¬ 
domly sampled from a uniform distribution between 0 
and 27r as to remove any azimuthal dependence at Tenc- 
Modelling the satellite system as an ensemble of test par¬ 
ticles ensure that the orbital evolution of the satellites 
is governed solely by gravitational interactions with the 
host planet and perturbations from the IG and the Sun. 

The semimajor axes of the satellites Os are initialized to 
±1% of the current orbital radius of Gallisto or lapetus. 
This fractional deviation is chosen to be on the order of 
the observed satellite eccentricities Cg (see Table 1). Be¬ 
cause changes in Og and Cg are related through changes 
in the satellites’ angular momentum, it is unlikely that 
larger shifts in Og would be reconcilable with the observed 
Cg assuming that each satellite formed out of a circum- 
planetary disk on a circular orbit (Ganup & Ward 2002; 
Mosqueira & Estrada 2003). Satellites are initialized on 
nearly coplanar and circular orbits {i < 0.1°,e < 10“®) 
as expected from circumplanetary formation scenarios. 
Satellite particles are not allocated physical sizes as we 
do not account for particle collisions in our simulations. 
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Fig. 2.— Distributions of the IG phase angle (0; bottom-left) and 
relative planet velocity (top-right) at tenc for encounters between 
Jupiter and an unspecified IG that resulted in the ejection of the 
latter. Relative planet velocities are uniquely determined by the 
impact parameter h and the specified fraction of the escape speed 
from Jupiter. Bottom-right: 2D histogram exhibiting the correla¬ 
tion between phase angle and relative velocity. The darkness of 
each bin is indicative of the logarithmic number of successful ejec¬ 
tions. The vertical dotted line is indicative of the escape velocity 
from the Sun at the semimajor axis of Jimiter. 

3.3. Numerical Code 

We performed our simulations using the REBOUND N- 
body numerical code (Rein & Liu 2012). We employ 
the IAS15 integrator (Integrator and Adaptive Step-size 
control, 15th order; Rein & Spiegel 2015) whose adap¬ 
tive timestepping ensures optimal resolution of the short 
satellite orbital periods and rapid encounter timescales.® 

4. RESULTS OF ICE GIANT EJECTIONS BY JUPITER 

Here we focus on simulations in which Jupiter is the 
satellite-hosting gas giant planet that ejects the hypo¬ 
thetical fifth giant planet from the solar system. 

4.1. Properties of Encounters 

Fig. 2 summarizes the properties of the planetary en¬ 
counters which result in the ejection of the IG planet. We 
find Nsim = 278 such encounters. The greatest impact 
parameter we find capable of ejecting the IG is ~ 0.05 
AU. 

All successful ejections involve closest approach in the 
lower hemisphere of the xy-plane in Jupiter’s reference 
frame (180° < 9{tenc) < 360°; see Fig. 1). As familiar 
from spacecraft gravity assists, an IG trailing Jupiter at 
tenc will receive a positive velocity kick via their interac¬ 
tion. At these phase angles, the relative velocity vector is 
rotated by the encounter such that the IG’s inertial speed 
accelerates. Depending on b and Vrei, the encounter can 
potentially boost the IG to escape velocity from the so¬ 
lar system. From simple vector diagrams, the maximum 
increase to the IG velocity can be shown to occur when 
d{tenc) = 270°. Fig. 2 highlights this as the distribution 
of 6{tenc) peaks at 270° where the encounter geometry 
is most conducive to ejecting the IG. 

Low Vrei trajectories lead more often to ejection since 
Jupiter can more effectively deflect the IG’s path. We 

® A short video from REBOUND depicting a close planetary en¬ 
counter with a ring of regular satellites can be found here. 



Fig. 3.— Evolution of five sample Jovian satellites as Jupiter 
undergoes a close encounter with an unspecified IG. The encounter 
parameters at tenc'- impact parameter 6, relative planet velocity 
Vnei, and IG phase angle 9 are shown in the legend of the lower 
panel. Each satellite’s azimuthal position relative to d{tcnc) is 
annotated in the upper panel. Vertical and horizontal dotted lines 
indicate the epoch of encounter and the current values of Callisto’s 
orbital elements respectively. 

note that there are many encounters in which the veloc¬ 
ity of the IG with respect to Jupiter at tenc is greater 
than the escape speed from the solar system at Jupiter’s 
distance from the Sun. This is due to the IG being ac¬ 
celerated upon approach to the encounter which occurs 
deep in Jupiter’s gravitational well. Such cases consist of 
the IG becoming unbound from the solar system prior to 
the realization of the impact parameter at tenc and are 
thus accelerated to super-escape speeds even before its 
closest approach to Jupiter. 

The 2D histogram in Fig. 2 depicts the correlation 
between the relative velocity of the planets at the epoch 
of encounter and the phase angle in the jovian-centric 
reference frame. Encounters which occur with the IG 
at a smaller heliocentric distance than Jupiter [Oitenc) < 
270°) require larger relative velocities in order to eject the 
IG. As the IG moves outwards (increasing 9) the relative 
velocities necessary to eject the IG decrease because the 
pull applied by Jupiter following the gravity assist is less 
efficient at decelerating the IG to sub-escape speeds given 
the IG’s trajectory. 

4.2. Resulting Jovian Satellite Orbits 

Glose planetary encounters between Jupiter and the 
IG can result in significant perturbations to the orbits 
of the in-situ Jovian satellites. The time evolution of Os 
and Cs for five sample Jovian satellites throughout one 
encounter simulation demonstrates this fact (Fig. 3). 

The five sample satellites which are initialized on 
nearly circular orbits get kicked to moderately eccentric 
orbits by a planetary encounter with b ~ 0.042 AU and 
Vrei ~ 13.57 km/s. Hence, the characteristic encounter 
timescale is b/vrei 5.3 days or about one third of Gal- 
listo’s orbital period. The orbital perturbation therefore 
acts roughly as an impulse to the satellites making their 
azimuthal position at tenc & parameter of importance. 
Indeed, in Fig. 3 the satellites whose azimuthal posi¬ 
tions are increasingly distant from 9(tenc) stre on average 
less perturbed than satellites which are close to 9{tenc)- 
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Semimajor axis (Rj) 

Fig. 4.— 2D histogram of the final average semimajor axes and 
eccentricities of Jovian satellites. The colorbar is indicative of the 
logarithmic number of satellites in each bin. The current values 
of Callisto’s semimajor axis and eccentricity are depicted on the 
plot as a white diamond. The horizontal dashed line highlights 
2ecallisto which marks the boundary between satellites that can 
(es < 2ecaiiisto) and cannot (e^ > 2ecallisto) be reconciled with 
Callisto’s current orbit. 

In addition, the position of the satellite relative to 0{tenc) 
determines the direction of the satellite’s radial shift. For 
example, satellites immediately trailing the IG at tenc 
will be pulled forward in their orbits (e.g. red curve) in¬ 
creasing Og, whereas satellites immediately upstream of 
the IG at tenc get pulled downward, thus decreasing Ug 
(e.g. black curve). 

The average final orbital elements of all Jovian satel¬ 
lites are depicted in Fig. 4 and show close agreement 
with the results from Deienno et al. (2014) despite the 
broader investigation of the close encounter parameter 
space. Ghanges in Ug and Cg are coupled through the 
satellite’s orbital angular momentum. This effect is high¬ 
lighted in Fig. 4, as satellites which are radially trans¬ 
ferred to either a significantly larger or smaller semimajor 
axis are those which exhibit the largest deviation from a 
circular orbit. 

For each satellite we define a reconcilable orbit to 
be when the satellite’s final average eccentricity is < 
‘^ecaiiisto (recall ecaiHsto ~ 0.007). The factor of 2 in our 
definition comes from the subsequent eccentricity evo¬ 
lution due to tidal damping of Gallisto (Deienno et al. 
2014, c.f. Fig. 7) and allows for Gallisto to be excited 
beyond the present ecaiUsto at tenc and consequently set¬ 
tle into its current orbital eccentricity in the 4 Gyrs fol¬ 
lowing the solar system’s instability phase. 

For the remainder of the paper we refer to the event of 
a simulated satellite’s final orbit being reconcilable with 
Gallisto as RS for “reconcilable satellite”. The boundary 
dividing RS from non-i?^ satellites is depicted in Fig. 4 
as a dashed horizontal line. 

Final average Cg values are never > 0.4, implying that 
no Jovian satellite becomes unbound from Jupiter follow¬ 
ing the ejection of the IG. That is, the vast majority of 
planetary encounters that are capable of ejecting an IG 
from the solar system are not sufficiently violent to strip 
Jupiter’s regular satellites. This favours the possibility 
of Jupiter being able to eject an IG whilst retaining a 
regular satellite whose orbit is Gallisto-like. We estimate 


the likelihood in Sect. 6. 

While the phase angle 9{tenc) is an important param¬ 
eter in determining whether the IG is ejected, it has lit¬ 
tle effect on the fraction of perturbed satellites due to 
their uniformly distributed azimuthal positions around 
Jupiter. Therefore, for each b and Vrei we can marginalize 
over 9 to show in Fig. 5 the fraction of satellites reconcil¬ 
able with the orbit of Gallisto after each encounter. Inter¬ 
polating over the fraction of reconcilable Jovian satellites 
in {b,Vrei) space, the resulting high resolution contours at 
10% and 50% are fitted with cubic functions and over¬ 
plotted as solid and dashed curves respectively to aid 
in visualization of the different regions of the parameter 
space. 

It is clear that as encounters become closer, the per¬ 
turbation to eg increases and the fraction of reconcilable 
satellites shrinks. Similarly, as the duration of encounters 
becomes longer (smaller Vrei), the timescale over which 
the perturbation is applied grows, thus increasing the de¬ 
viation of satellite orbits from circular. Hence, forming 
and maintaining a Jovian satellite with a Gallisto-like 
orbit favours encounters which are wide and fast. 

We suggest that researchers simulating solar system 
formation scenarios can use Fig. 5 to estimate whether 
or not a given Jupiter/IG encounter is consistent with 
Jupiter’s current Galilean satellite orbits. We note that 
our requirement that the IG gets ejected only determines 
what regions of the plot are populated, hence the con¬ 
tours are useful guides regardless of whether or not the 
IG survives the encounter (and for any 9). We caution, 
however, that these will only be approximate due to our 
assumption of coplanarity. For a given b and Vrei, in¬ 
clined encounters will increase the fraction of satellites 
whose eccentricities are reconcilable with current orbits. 
However, they will also tend to overly excite the satel¬ 
lite inclinations. We discuss this further in Sect. 7.1. 
For particular simulations that lie on the boundary of 
plausibility in Fig. 5, one could re-simulate the close 
approaches, including inclined encounters, following our 
setup to more accurately quantify the effect. 

5. RESULTS OF ICE GIANT EJECTIONS BY SATURN 

Here we focus on simulations in which Saturn is the 
satellite-hosting gas giant planet that ejects the hypo¬ 
thetical fifth giant planet from the solar system. 

5.1. Properties of Encounters 

The properties of planetary encounters between Saturn 
and an IG are summarized in Fig. 6. We find Ngim = 274 
simulations which result in the IG being ejected by Sat¬ 
urn. The largest impact parameter that is still capable 
of ejecting the IG is nearly identical to the Jupiter/IG 
case; ss 0.05 AU. We perform an additional sampling of 
encounter parameters with h > 0.05 AU at an increased 
resolution to confirm that the largest impact parameter 
capable of ejecting the IG is nearly the same in both the 
Jupiter and Saturn cases rather than being a statistical 
anomaly due to the finite sampling of encounter param¬ 
eters. We find that this limit is real and not an artifact 
of our sampling procedure. The physical interpretation 
of this similarity is beyond the scope of this paper. 

Similar trends to those shown in Fig. 2 are observed. 
This should be expected because, modulo the varia¬ 
tions in the physical parameters Mp and a^, the Sat- 
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Fig. 5. — The fraction of reconcilable Jovian satellites as a function of the encounter’s impact parameter and relative planet velocity. 
The colorbar indicates the fraction of satellites whose final average eccentricity is less than or equal to ^ecalHsto- Smoothed contours at 
10% (solid line) and 50% (dashed line) assist in visualizing the different regions of the parameter space. The survival fractions shown here 
reveal that Jupiter has a reasonably large probability of ejecting an IG while retaining a Callisto-like satellite (see Sect. 6). 



Fig. 6.— Same as Fig. 2 for encounters between Saturn and an 
unspecified IG. The vertical dotted line is indicative of the escape 
velocity from the Sun at the semimajor axis of Saturn. 

urn/IG encounters investigated are fundamentally equiv¬ 
alent to those described in Sect. 4.1. With the excep¬ 
tion of the decreased magnitude of Vrei by a factor of 
\/Msat /Mjup on average, the set of successful scattered- 
by-Saturn simulations are statistically identical to those 
in the scattered-by-Jupiter simulations; i.e. the his¬ 
tograms in Figs. 2 and 6 exhibit the same behaviour. 

5.2. Resulting Kronian Satellite Orbits 

Similarly to the case of Jupiter ejecting the IG from 
the solar system, the Kronian satellite orbits will be per¬ 
turbed as a result of the encounter between Saturn and 
the IG. As a result of the wide separation of lapetus 
61 Rs) and the similar nature of the encounters be¬ 


tween an IG and either Jupiter or Saturn (Figs. 2 & 6), it 
is expected that such orbital perturbations will be more 
destructive to simulated satellite orbits than to those pre¬ 
viously explored in Sect. 4. This notion is supported in 
Figs. 7, 8, and 9 when compared to their counterparts in 
the Jupiter/IG case (Figs. 3, 4, and 5) as the fraction 
of reconcilable Kronian satellites is systematically lower 
than in the Jupiter case. 

The time evolution of Og and Cg for five example satel¬ 
lites during a sample planetary encounter is shown in 
Fig. 7. This example is a particularly violent encounter 
with b ~ 0.033 AU and Vrei ~ 5-87 km/s as no satellite in 
this simulation has a resulting orbit that is reconcilable 
with lapetus. One satellite, whose azimuthal position 
is only 4.6° from 0{tenc), is ejected from Saturn by the 
encounter. The encounter timescale {bfvrei 9.6 days) 
is much less than the orbital period of lapetus making 
the effect of the encounter behave approximately as an 
impulse. 

The average final orbital elements of all Kronian satel¬ 
lites after each unique encounter are depicted in Fig. 
8. The black contours in Fig. 8 are representative of 
the final orbital elements of the Jovian satellites from 
Fig. 4 to aid in visual comparison between the final 
orbits of Jovian satellites to Kronian satellites follow¬ 
ing similar close planetary encounters. This shows that 
the wide-separation of the simulated Kronian satellites 
makes them more susceptible to excessive orbital pertur¬ 
bations. 

The event RS in the Saturn/IG case is defined similarly 
to the Jupiter/IG case and is achieved by Kronian satel¬ 
lites whose final average Cg < eiapetus- As noted earlier, 
the eccentricity evolution of lapetus since the solar sys¬ 
tem’s instability phase is negligible (Gastillo-Rogez et al. 
2007) thus we assume that resulting orbits from our sim- 
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Fig. 7.— Same as Fig. 3 featuring the evolution of five sample 
Kronian satellites as Saturn undergoes a close encounter with an 
unspecified IG. Vertical and horizontal dotted lines indicate the 
epoch of encounter and the current values of lapetus’ orbital ele¬ 
ments respectively. 




Fig. 8. — Similar to Fig. 4 for the Kronian satellites with the in¬ 
clusion of Saturn’s Hill radius (~ 1100 Rs \ vertical dotted line) and 
contours of the 2D histogram shown for Jupiter (Fig. 4) for direct 
comparison. The horizontal dashed line highlights ejape.tua which 
marks the boundary between satellites that can (es < eiapetua) 
and cannot {cs > eiapetus) be reconciled with lapetus’ current 
orbit. Note the logarithmic x-axis. 


ulations will go largely unchanged to the present day. 
The boundary dividing RS from non-i?S' is depicted in 
Fig. 8 as a dashed horizontal line. 

As noted in Sect. 4.2, no Jovian satellite becomes more 
eccentric than 0.4 whereas a sizable fraction of Kronian 
satellites become equivalently or excessively excited in¬ 
cluding a subset of Kronian satellites whose final av¬ 
erage Cs is approximately unity. Furthermore, we find 
that ~ 6% of sampled Kronian satellites are sufficiently 
excited to final Cs > 1. The orbital elements of these 
ejected satellites are not included in Fig. 8. It is clear 
that the perturbations to the Kronian satellites resulting 
from close encounters, are on average much stronger and 
capable of stripping lapetus-like satellites from the Kro¬ 
nian system. The large fraction of satellites perturbed 
beyond the orbit of lapetus makes it statistically diffi¬ 


cult for Saturn to have ejected an IG whilst retaining an 
lapetus-like regular satellite. We formally estimate the 
likelihood of reconciling lapetus’ orbit in Sect. 6. 

The effect of encounter properties on the resulting 
Kronian satellite orbits is shown in Fig. 9. For each 
encounter with a given impact parameter and relative 
planet velocity, the fraction of Kronian satellites whose 
final orbit is reconcilable with the current orbit of lapetus 
is shown. Approximate contours of 10% and 50% frac¬ 
tions are over-plotted. The 10% contour is computed 
identically to the contours in Fig. 5 (i.e. cubic inter¬ 
polation). However, the 50% contour lacks a sufficient 
number of points to perform a robust cubic interpola¬ 
tion. Therefore we opt for a linear fit in its place. It is 
clear that the region of the {b,Vrei) parameter space in 
which > 50% of the Kronian satellites are reconcilable 
with lapetus’ orbit is very small compared to the Jovian 
satellites with only five unique ejections sampled there. 

The effects of b and Vrei on resulting satellite orbits are 
nearly identical to those shown in Fig. 5. One interest¬ 
ing difference which is unique to Saturn/IG encounters 
is that regardless of how close the encounter is, if the en¬ 
counter is sufficiently long {vrei ^ 8 km/s), then < 10% 
of satellites will be reconcilable with lapetus. This is in 
contrast to the Jovian satellite case where even the slow¬ 
est encounters could retain a high fraction of Gallisto-like 
satellites if the encounter’s impact parameter is large. 
Another important feature to note is that no Saturn/IG 
encounter resulting in the ejection of the latter is guar¬ 
anteed to preserve an lapetus-like satellite. That is, even 
the least violent encounters leading to ejection are only 
capable of preserving a maximum of 70% of the in-situ 
lapetus-like satellites. 

Fig. 9 contains information on the likelihood that la¬ 
petus survives an IG ejection by Saturn with a particu¬ 
lar b{tenc) and Vrei{tenc) when the mutual inclination of 
the encounter is near zero. Similarly to the Jupiter/IG 
encounter case, we suggest that researchers simulating 
solar system formation scenarios can use Fig. 9 to es¬ 
timate whether or not a given Saturn/IG encounter is 
consistent with lapetus’ current orbit in the limit of un¬ 
inclined planetary encounters. 

6. LIKELIHOOD OF RECONCILING SATELLITE ORBITS 
FOLLOWING ICE GIANT EJECTIONS 

6.1. Methodology 

The resemblance of the final average orbits of the Jo¬ 
vian and Kronian satellites (Figs. 4 and 8) to the current 
orbits of Callisto and lapetus can, in principle, be used 
to compute the likelihood of a fifth giant planet getting 
ejected by either of the gas giant planets in the early so¬ 
lar system. If such an event were to have occurred, it 
must be consistent with the satellite orbits presently ob¬ 
served. The likelihood of an IG getting ejected by either 
gas giant requires the likelihood of the current orbits of 
Gallisto or lapetus being reconciled by simulated satel¬ 
lites after an IG ejection. A successful event in which the 
resulting orbit of a simulated Jovian (Kronian) satellite 
is reconcilable with Callisto’s (lapetus’) current orbit is 
referred to as RS for “reconcilable satellite”. 

Given that an IG gets ejected in all simulations (event 
IGE] “ice giant ejected”), for the satellite in the 
simulation, we record whether or not the event RS is 
achieved by 
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Fig. 9.— Same as Fig. 5 for the Kronian satellites. The colorbar indicates the fraction of satellites whose final average eg is less than or 
equal to ejapetus- The vertical dotted line indicates lapetus’s current semimajor axis. The survival fractions here reveal that Saturn has a 
small probability of ejecting an IG while retaining an lapetus-like satellite if we assume that lapetus formed like a regular satellite out of 
a circumplanetary disk. 


p.,(fiS|/G£;) = { J (1) 

However, not all planet orbits in our ejection simula¬ 
tions are statistically relevant. Our adopted methodol¬ 
ogy for determining which encounter parameters result 
in an ejection is heavily biased towards initially high- 
eccentricity orbits of the IG (e/ < 1) which are not long¬ 
term stable and therefore are uncommon in nature. This 
bias naturally arises because it is easier to kick the IG to 
e/ > I if 6/ is initially very close to unity. Therefore in 
each simulation j, each satellite’s likelihood of RS must 
be weighted by the likelihood of the IG having an ini¬ 
tially bound orbit with e/j, where eij is the IG’s initial 
eccentricity in the j*** simulation. 

The corresponding weighting function W{eij) is mod¬ 
elled by the distribution of planet eccentricities observed 
in a statistically significant number of exoplanetary sys¬ 
tems as the eccentricity distribution of solar system bod¬ 
ies is insufficient as it can only be derived in the limit 
of small number statistics. These data are recovered 
from the www.exoplanets.org database (Han et al. 2014) 
and only include RV exoplanet detections with Doppler 
variation semi-amplitude K /> 5 (neglect low signal- 
to-noise observations). We find that our results do not 
sensitively depend on whether we focus on RV detec¬ 
tions or the full catalogue of exoplanets with orbital so¬ 
lutions. The distribution of planet eccentricities is shown 
in Fig. 10. Due to variations in empirically derived ec¬ 
centricity distributions we consider three proposed ana¬ 
lytical forms of the distribution. Namely a Beta proba¬ 
bility density function (PDF) (Kipping 2013), a Rayleigh 
PDF plus decaying exponential (Juric & Tremaine 2008; 
Steffen et al. 2010), and the model from Shen & Turner 



Fig. 10. — Histogram of the planet eccentricities used to model 
the weighting function W{ej k) in Eqs. 2 and 3. The three analyt¬ 
ical models used to fit the distribution are overplotted (solid line: 
Beta PDF, dotted line: Rayleigh PDF plus exponential, dashed 
line: model from Shen Sz Turner (2008); see text). 

(2008). We use a Levenberg-Marquardt least squares 
algorithm to fit for each PDF’s unique parameters. A 
summary of the adopted distributions including fitted 
parameters is given in Table 2 and each fitted weighting 
function is over-plotted in Fig. 10. Differences in W{eij) 
from adopting three unique PDFs results in ^ 6% vari¬ 
ance among computed likelihoods. 

While we take the known exoplanet orbital eccentrici¬ 
ties as our nominal distribution, this sample is obviously 
biased by the requirement that systems be stable. If 
early on, the solar systems ice giants moved on signifi¬ 
cantly more elliptical paths than the known exoplanets, it 
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TABLE 2 



Analytical fits to the planet eccentricity 

DISTRIBUTION. 



PDF Name 

PDF 

Parameter 1 

Parameter 2 

Parameter 3 

Beta 

Rayleigh + exponential 

ST08 

Wfi{x-,a,b) %1 x)^ ^ 

ITR+exp(a;;a, A,(t) = exp {^) + aXe^p {-\x) 

WsT 0 s(x-,k,a)- 2 ") 

0.786 ± 0.055 

0.782 ± 0.301 

0.305 ± 0.012 

2.764 ±0.167 

5.131 ±2.342 

3.413 ± 0.309 

0.266 ±0.061 


Note. — Parameter columns are written in the order in which they appear in the functional form W{x -,..shown in the PDF column. 
We write the independent variable of planet eccentricity as x to make the distinction between eccentricity and Euler’s number e ^ 2.7182818. 


would be easier to eject planets while keeping Jupiter and 
Saturns satellites on nearly circular orbits, raising the 
likelihood of reconciling satellite orbits. However, this 
scenario would require a substantial amount of damping 
to subsequently recircularize the surviving planets orbits, 
and if eccentricities reached such high values, one might 
expect to also lose Uranus and Neptune. These consid¬ 
erations may be interesting directions for future work. 

One must additionally weight pij by the simulation’s 
impact parameter bj , since wider encounters should occur 
more frequently. Given our coplanar setup, the differen¬ 
tial interaction cross-section for a given impact parame¬ 
ter scales linearly with bj rather than with b^ as is true 
in the full 3D case. We find that our results do not sen¬ 
sitively depend on this distinction. We therefore assume 
W'{bj) oc ^ 

The distribution of ejj is not explicitly prescribed 
and is instead determined from the encounter parameters 
(b,Vrei-,&)- The resulting distribution of ejj is not sam¬ 
pled uniformly, unlike the distribution in 6 or 0, implying 
that the step size Ae/ when integrating over IG eccen¬ 
tricities is not constant throughout the domain. Hence, 
simulations sample IG eccentricity bins of various widths 
which are taken into account when computing the likeli¬ 
hood function by calculating the width of IG eccentric¬ 
ity bins among the Wim simulations. In this way, the 
subspace of IG eccentricities which is over-sampled gets 
averaged over. Similar factors Abj and A0j are included 
but are constants, because the parameters are sampled 
uniformly, and therefore do not affect the resulting like¬ 
lihood. 

Gombining the aforementioned effects into the likeli¬ 
hood of obtaining a reconcilable satellite {RS) orbit given 
an IG ejection (IGE), we write 

P{RS\IGE) = 

, N Nsim 

-EE P^ARS\IGE)WieI,,)W'ib,)Aep,Ab,Ae, 

^ i=i j=i 

( 2 ) 

where 

Nairn 

fi = NY^ W{epj)W'{b,)Aei,,Ab,Ae„ (3) 
i=i 

is the normalization factor of the weighted mean and 
we are careful to account for the fact that our uniform 
sampling of {b,Vrei,0) leads to a non-uniform distribution 
of eij. 

6.2. Likelihoods 


From our Nsim = 278 ejection simulations by Jupiter, 
each with N satellites, we calculate that the likelihood 
of Jovian satellite orbits remaining consistent with the 
observed orbit of Callisto, is 42%. This value is the 
median of the results we obtained from adopting the 
three distinct eccentricity weighting functions discussed 
in Sect. 6.1 (see Table 2). The median absolute deviation 
among the likelihoods is small; < 1%. 

By contrast, in our Ngim = 274 ejection simulations by 
Saturn, the likelihood of Kronian satellite orbits remain¬ 
ing consistent with lapetus’ observed orbit is ^ 1%, more 
than an order of magnitude less likely than for Gallisto 
around Jupiter. The main reason for this wide likelihood 
disparity is that a given encounter will perturb lapetus 
more strongly than Callisto because the former is less 
tightly bound to Saturn than the latter is to Jupiter. 
In this case, the median absolute deviation among the 
three eccentricity weighting functions is < 0.6% which is 
comparable to the likelihood itself. 

7. DISCUSSION 

To recapitulate, the above likelihoods assumed an IG 
ejection in a coplanar geometry, considering whether an 
initially circular Callisto (lapetus) would have its orbital 
eccentricity excited beyond values reconcilable with its 
current orbit. We now consider these assumptions in 
turn, in order to interpret the results from our study and 
discuss possible implications of our work. 

7.1. Effect of Inclined Encounters 

To limit the computational cost of our study, we have 
restricted our analysis to planetary encounters with no 
mutual inclination. However, if planetary eccentricities 
are excited enough to permit orbit-crossing and ejections, 
one might expect comparably large orbital inclinations. 
During an inclined encounter, some of the applied torque 
goes into realigning Gallisto’s (lapetus’) orbital plane so 
that, on average, the perturbations to Os and Cg are re¬ 
duced as some energy goes into increasing ig. 

A preliminary analysis of inclined encounters with 
mutual inclination iftenc) = 5° effectively revealed no 
change to the final satellite orbital elements compared 
to uninclined, but otherwise equivalent, encounters. In 
a more heavily inclined test case with i{tenc) = 45°, 
we found that, on average, the final Cg were smaller by 
> 0.01 than in the uninclined case. 

Therefore, by limiting our investigation to coplanar 
encounters, we are sampling the largest possible per¬ 
turbations to Cg without affecting ig. The introduction 
of inclined encounters would thus raise the likelihoods 
quoted in Sect. 6.2. However, such a 3D case would ad¬ 
ditionally excite the satellite inclinations. As mentioned 
by Deienno et al. (2014), these inclinations may provide 
more rigorous constraints on planetary encounters be- 
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cause even in cases where eccentricity damping is impor¬ 
tant, the tidal evolution of the inclinations is effectively 
null. A generalization of this study would therefore con¬ 
sider the combined constraint from the satellites’ orbital 
eccentricities and inclinations. However, the uncertain¬ 
ties in the likelihoods we compute in Sect. 6.2 are domi¬ 
nated by the large uncertainties in the planetary orbital 
eccentricities (and inclinations) early in the solar system. 
We therefore believe that our simplified analysis consid¬ 
ering only the orbital eccentricities captures the correct 
likelihoods at the approximate level allowable by our cur¬ 
rent state of knowledge. 

As a rough check, in a fully 3D case, the differential 
cross-section for encounters of a given impact parameter 
would scale as rather than as bj. Using the full data 
set from our coplanar simulations, but adopting this new 
scaling for W'{bj) in Eqs. 2 and 3, we find that our results 
do not change discernibly. Specifically, in the Jupiter 
case our results change from ^ 42% to 54% and from ^ 
1.1% to 1.3% in the Saturn case. This is certainly within 
the errors of our uncertain knowledge of the planetary 
orbits’ initial eccentricities and inclinations and therefore 
does not appreciably change the interpretation of our 
results. 

7.2. Interpretation of Results 

We conclude that Jupiter could plausibly have ejected 
an IG from the solar system. Nearly half of the hypothet¬ 
ical set of ejections that were modelled in Sect. 6.2 keep 
Callisto on an orbit which is reconcilable with the one 
we observe today. Put another way, we conclude that 
Callisto’s orbit cannot meaningfully constrain whether 
or not Jupiter ejected an additional IG in the early So¬ 
lar System. Nevertheless, this is an important test for 
the fifth-giant-planet hypothesis to pass, thus providing 
more stringent evidence for its plausibility. 

Interpretation of our results in the Saturn case is more 
subtle. We showed in Sect. 6.2 that an initially circu¬ 
lar lapetus orbit gets overly excited by a single ejection 
event ^ 99% of the time. This suggests that Saturn is 
not capable of ejecting an IG mass planet from the solar 
system. 

But did lapetus originally move on a circular path, 
as expected if it formed out of a circumplanetary disk? 
Starting with a circular orbit, all encounters act to raise 
the eccentricity. But if the moon’s initial path were in¬ 
stead elliptical, some ejection geometries could act to 
lower the eccentricity, complicating the constraint that 
we nominally set. Perhaps one reason to doubt that la¬ 
petus formed from a circum-kronian disk is that it is sub¬ 
stantially inclined (« 8°) to the local equilibrium plane 
that one would expect it to follow. 

So what caused this aberrant inclination? Hamilton 
(2013) recently suggested that lapetus could indeed have 
formed on a circular orbit from a disk in its equilibrium 
plane, and subsequent collisions between Saturn’s inner 
moons could have instead tilted the equilibrium plane for 
the exterior moons by the requisite amount. In this case, 
our study implies that Saturn likely did not eject an IG 
from the solar system because lapetus’ relatively low ec¬ 
centricity orbit cannot be reconciled with IG ejections by 
Saturn. On the other hand, lapetus’ inclination could be 
the signature of an ejection event itself. Nesvorny et al. 
(2014) recently studied such a scenario trying to simulta¬ 


neously match lapetus’ current orbital eccentricity and 
inclination using simulations of the specific early solar 
system dynamical instabilities from NM12. They show 
that some cases are capable of sufficiently exciting lape¬ 
tus’ orbital inclination whilst maintaining a low orbital 
eccentricity even for encounters as close as those consid¬ 
ered in this study but not necessarily leading to ejection 
of the IG. 

In the Kronian case, therefore, the interpretation de¬ 
pends critically on the formation mechanism for lapetus, 
which is currently unknown. If lapetus’ inclination is 
the result of collisions between inner moons (Hamilton 
2013), our results show that Saturn is unlikely to have 
ejected an IG from the early Solar System. By contrast, 
if lapetus’ orbit is fully explainable through close plane¬ 
tary encounters, Nesvorny et al. (2014) showed that this 
requires many such close approaches which might be ca¬ 
pable of ejecting the IG. 

7.2.1. Single-Encounter Assumption 

We argue that the sole consideration of the final gas 
giant/lG encounter leading to ejection of the latter, is a 
sufficient measure of how close planetary encounters will 
modify regular satellite orbits. Nesvorny et al. (2014) 
showed for lapetus, which is more susceptible to dynam¬ 
ical perturbations than Callisto, that numerous ‘soft’ en¬ 
counters prior to ejection have a fractional effect on satel¬ 
lite eccentricity compared to the eccentricity kick typi¬ 
cally felt during the final encounter (Fig. 4 and 8). It 
should be noted that inclusion of numerous ‘soft’ encoun¬ 
ters prior to ejection would be detrimental to the com¬ 
puted likelihoods P{RS\IGE) thus emphasizing that our 
results represent a conservative, best-case scenario. 

7.2.2. Close Encounters Without Ejection 

Due to the phase angle of an encounter being ran¬ 
dom, it is possible for a close encounter to have oc¬ 
curred at 0° < 0 < 180° and is subsequently not in¬ 
cluded in our likelihood calculation as all ejection events 
occurred at 180° < 9 < 360° (see Figs. 2 and 6). Any 
close encounter that does not lead to an ejection would 
still perturb regular satellite orbits, making their orbits 
prior to an IG ejection event, non-circular. The effect 
of an increased initial eccentricity reduces P[RS\IGE), 
again making our calculation a conservative upper limit. 
However, because the frequency of close encounters with 
0° < 0 < 180° is equivalent to the frequency of close en¬ 
counters with 180° < 9 < 360°, the effect on P{RS\IGE) 
of close encounters without an ejection will only differ 
from our calculated values by a factor of order unity, 
such that we capture the correct order of magnitude on 
P{RS\IGE). 

7.3. Additional Considerations 

One sought after quantity relating to this work is the 
probability of an IG getting ejected by either gas giant 
given that the orbit of Gallisto or lapetus can be rec¬ 
onciled; P{IGE\RS). Using Bayes’ theorem, this quan¬ 
tity could in principle be computed with knowledge of 
the likelihood functions calculated in this paper but also 
requires the independent probability of the IG getting 
ejected {P{IGE); see NM12) and the normalization fac¬ 
tor by the probability of a satellite exhibiting the current 
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orbit of one of the wide-separation satellites {P{RS)). 
Despite the abundance of work done on the latter 
(e.g. Canup & Ward 2002; Estrada & Mosqueira 2006; 
Ward & Canup 2010; Crida & Charnoz 2012; Hamilton 
2013; Heller et al. 2015), precise constraints on P{RS) 
are difficult to compute. 

Also, a robust calculation of the likelihood of ejecting 
an IG from the solar system would require the event to 
be consistent with a vast number of constraints imposed 
by solar system bodies of which the orbital constraint 
RS imposed by Callisto or lapetus is just one. There¬ 
fore, such a calculation is not practical. However, con¬ 
sistency checks of IG ejections by a gas giant, such as 
those presented in this study, help to substantiate the 
proposed existence of a fifth giant planet. Based on our 
current understanding, which includes the main results 
of our study, there is little evidence demonstrating that 
an additional IG mass planet could not have existed in 
the early solar system. 

8. SUMMARY 

Several studies trying to match the solar system’s cur¬ 
rent orbital architecture argue for an early period of 
frequent planetary encounters (e.g. Tsiganis et al. 2005; 
Brasser et al. 2009; Morbidelli et al. 2009; Levison et al. 
2011). In addition, Nesvorny (2011) found that adding 
a fifth giant planet to the solar system, which is sub¬ 
sequently ejected, better matches the current orbits of 
the remaining giant planets. In this paper, we therefore 
study whether such an ejection by either Jupiter or Sat¬ 
urn is reconcilable with the current observed orbits of 
their outermost regular satellites, Callisto and lapetus. 
Our main conclusions are as follows: 

• The properties of planetary encounters (i.e. impact 
parameter, relative planet velocity, and encounter 


Batygin, K., Brown, M. E., & Betts, H. 2012, ApJ, 744, L3 
Brasser, R., Morbidelli, A., Gomes, R., Tsiganis, K., &; Levison, 
H. F. 2009, A&A, 507, 1053 
Canup, R. M., & Ward, W. R. 2002, AJ, 124, 3404 
Castillo-Rogez, J. C., Matson, D. L., Sotin, C., et al. 2007, Icarus, 
190, 179 

Crida, A., & Charnoz, S. 2012, Science, 338, 1196 
Deienno, R., Nesvorny, D., Vokrouhlicky, D., & Yokoyama, T. 
2014, AJ, 148, 25 

Delorme, P., Gagne, J., Malo, L., et al. 2012, A&A, 548, A26 
Estrada, P. R., & Mosqueira, I. 2006, Icarus, 181, 486 
Fernandez, J. A., & Ip, W.-H. 1984, Icarus, 58, 109 
Gomes, R., Levison, H. F., Tsiganis, K., & Morbidelli, A. 2005, 
Nature, 435, 466 

Greenberg, R., & Van Laerhoven, C. 2011, ApJ, 733, 8 
Hahn, J. M., & Malhotra, R. 1999, AJ, 117, 3041 
Hamilton, D. P. 2013, in AAS/Division for Planetary Sciences 
Meeting Abstracts, Vol. 45, AAS/Division for Planetary 
Sciences Meeting Abstracts, 302.01 
Han, E., Wang, S. X., Wright, J. T., et al. 2014, PASP, 126, 827 
Heller, R., Marleau, G.-D., & Egon Pudritz, R. 2015, ArXiv 
e-prints, arXiv:1506.01024 
Juric, M., & Tremaine, S. 2008, ApJ, 686, 603 
Kipping, D. M. 2013, MNRAS, 434, L51 

Levison, H. F., Bottke, W. F., Gounelle, M., et al. 2009, Nature, 
460, 364 

Levison, H. F., Morbidelli, A., Tsiganis, K., Nesvorny, D., & 
Gomes, R. 2011, AJ, 142, 152 

Levison, H. F., Morbidelli, A., Van Laerhoven, C., Gomes, R., & 
Tsiganis, K. 2008, Icarus, 196, 258 


geometry) between Jupiter or Saturn and an un¬ 
specified ice giant, which are sufficiently violent to 
eject the latter, exhibit similar trends. 

• The current (dynamically cold) orbit of the widest- 
separation Galilean satellite Callisto, has a signifi¬ 
cant likelihood (^ 42%) of being reconciled follow¬ 
ing the ejection of an ice giant planet from the solar 
system by Jupiter. 

• Given the observed difficulty in reconciling the or¬ 
bit of lapetus with simulated Kronian satellites fol¬ 
lowing an ejection event, the likelihood of Saturn 
ejecting an ice giant from the solar system is deter¬ 
mined to be unlikely; likelihood ~ 1%. 

• However, we note the caveat that this interpreta¬ 
tion is heavily dependent on the assumed forma¬ 
tion scenario of lapetus out of circum-kronian disk. 
Currently, the formation of lapetus is largely un¬ 
certain. 

We caution that these likelihoods should not be inter¬ 
preted in an absolute sense. Rather, they are useful in 
showing that it is much easier for Jupiter to have ejected 
an ice giant than it is for Saturn. The evident plausibil¬ 
ity of Jupiter being able to eject an IG thus supports the 
hypothesis of a fifth giant planet’s existence in the early 
solar system. 
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